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Spatial and temporal evolution is studied of two powerful short laser pulses 

o . 

. having different wavelengths and interacting with a dense three-level Lambda- 

, type optical medium under coherent population trapping. A general case of 

O ■ 

, unequal oscillator strengths of the transitions is considered. Durations of the 

g: ................ <.>.,»e—...„. 

than any of the relevant atomic relaxation times. We propose analytical and 
numerical solutions of a self-consistent set of coupled Schrodinger equations 
and reduced wave equations in the adiabatic limit with the account of the first 



^ . 

^ ^, non-adiabatic correction. The adiabaticity criterion is also discussed with the 

account of the pulse propagation. The dynamics of propagation is found to 
be strongly dependent on the ratio of the transition oscillator strengths. It 
is shown that envelopes of the pulses slightly change throughout the medium 
length at the initial stage of propagation. This distance can be large compared 
to the one-photon resonant absorption length. Eventually, the probe pulse is 
completely reemitted into the coupling pulse during propagation. The effect 
of localization of the atomic coherence has been observed similar to the one 
predicted by Fleischhauer and Lukin (PRL, 84, 5094 (2000). 
PACS numbers: 42.50.Hz, 42.50.Rh, 42.65.Ky 
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I. INTRODUCTION 



The electromagnetically induced transparency (EIT) can be used to make optically thick 
media transparent to resonant laser radiation |]T[]. The EIT is the result of various quantum 
interference effects such as nonlinear interference coherent population trapping (CPT) 
adiabatic population transfer (APT) |^. Optical characteristics of the matter undergo 
drastic changes under those effects to such an extent that they can now be manipulated. A 
lot of interesting applications based on that have been proposed and experimentally realized 
(see, e.g. M^)- 

Interesting and unusual phenomena caused by the above indicated effects can be observed 
when laser pulses propagate in a resonant three-level medium. Propagation of pulses under 



EIT conditions was studied for example in |11-18|. As a rule, the situations are consid 



ered when both pulses have identical forms and their duration is longer than the relaxation 
time of the intermediate resonant state (matched pulses [|ll| ; dressed field pulses [0,|3 



or 



when the duration of the coupling radiation considerably exceeds that of the probe radiation 
(adiabatons [jl5|,|16|,|19|). A theoretical study of certain features of spatial evolution under 



APT conditions is presented in [116 i p!7 i P0i . Propagation of soliton-like pulses in a three-level 
system is studied in [^. A three- level system with equal oscillator strengths is under con- 
sideration in all the above mentioned studies, whereas in actual fact the transition oscillator 
strengths most often are different. 

In this paper, the spatial and temporal evolution is studied of two overlapping short 
laser pulses propagating in a resonant optically thick medium that consists of three-level 
A-atoms. Pulses of such configuration are widely used to enhance the efficiency of nonlinear 



generation processes p2|,^. The two pulses are assumed to have identical shapes but 
different durations (T2 > Ti) as shown in Fig.l. The pulse durations are much shorter than 
any of the times of relaxation in the medium (short pulses). It is also assumed that the 
pulse envelopes satisfy the adiabaticity criterion |T3| , p^ , |25[] : 

(G2G1 - GiG2)/G^ < 1, (1) 



where Gi^2 are the Rabi frequencies of the respective fields, G = yGf + G^] the dot refers to 
time derivatives. Condition (|l|) is easy to satisfy by making one of the pulse amplitudes or 
both of them large even for short pulses. This will induce a strong coherence at the Raman 
transition resulting in the effect of CPT. The latter considerably decreases the absorption 
of the propagating resonant pulses. The dynamics of propagation of such pulses is studied 
here without restriction of the relationship between oscillator strengths of transitions. 

Our theoretical model involves a set of coupled Schrodinger equations and a set of reduced 
wave equations, allowing thus a simultaneous description of temporal and spatial evolution 
of the atomic system and the radiation. The equations are analyzed in approximation (|l]) 
with the account of the first non-adiabatic correction. It will be shown that the dynamics of 
propagation strongly depends on the oscillator strengths ratio. Also the spatial and temporal 
behavior is analyzed of the atomic Raman coherence. A possibility of localizing the atomic 
coherence spatially has been established. The results obtained are compared to the results 



reported in 26 



The paper is organized as follows. In Section II, we describe the model and present 
the basic equations. Section III contains self-consistent solutions of those equations in the 
adiabatic limit and describes the temporal behaviour of the level populations and the atomic 
Raman coherence in the optically thin medium. The spatial evolution of pulses in the 
optically thick medium is described in Section IV for various oscillator strength ratios. In 
Section IV, we also discuss the adiabaticity criterion and demonstrate the effect of spatial 
localization of the atomic coherence. Finally we summarize the results obtained. 



II. BASIC EQUATIONS 

The three-level system under consideration is shown in Fig.l together with the temporal 
configuration of the pulses as they enter the medium. The pulses travel along the same di- 
rection, z. States |0), |1) and |2) are connected by laser pulses Ei = l/2Ei(t) exp[—i{uJit — 
kiz)] + c.c. and E2 = 1/2 E2{t) exp[—i{uj2t — k2z)] + c.c, respectively. In our further consid- 
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eration we shall refer to the first pulse Ei as the probe and to the second pulse E2 as the 
coupling pulse. The Rabi frequency of probe pulse is comparable with that of the coupling 
pulse. The pulses are sent simultaneously into an atom. The pulse durations Ti_2 (^2 > Ti) 
are assumed to be much less than any of the relaxation times of atoms. The transition 
|0) — |2) is electric dipole forbidden. The intermediate state |1) is in a one-photon resonance 
with each field interacting only with the respective transition. 

The following standard set of equations describes the spatial and temporal dynamics 
of the probability amplitudes of atomic states &o,i,2 and slowly varying Rabi frequencies 
Gi = dioEi{t)/2h, G2 = d2iE2{t) / 2h in the local-time coordinate system r = t — z/c: 

dbo 

dr 
dh2 



g iClbi exp {—ikiz) , 



iG*2biexp{-ik2z), 
dbi 

= iGibo exp {ikiz) + iG'2&2 exp {ik2z). (2) 
= zi^i6i6oexp {ikiz), 

oz 
dG2 



dz 



iK2bib*2 exp {ik2z) . (3) 



Here we assumed zero one-photon detunings. Ki^2 = ttcji {diQ^ufN/hc are the propa- 
gation coefficients, is the atomic concentration, (iio,i2 are the dipole transition matrix 
elements, 001^2, ^1,2 are the frequencies and the wave numbers of the interacting waves in 
vacuum, c is the light velocity in vacuum. All atoms are assumed to be initially in the 
ground state |0): bo{—oo,z) = 1, bi^2{—oo, z) = 0. We use Gaussian pulses at the medium 
entrance z = for the purpose of numerical simulation: G'i(r) = Cj* exp(— r^/n2/T^^), 
G2{t) = G'°exp[-r2/n2/T2]. 

In terms of = bo exp {ikiz), 02 = &2exp {ik2z), ai = ibi, equations (||) and (H) can be 
written as: 

dao _ 

— - G^ai, 



^ = -Giao - G2a2. (4) 
dGi 



dG2 



K2aial. (5) 



dz 

The coupled equations (^) and (^) give a complete semiclassical description of the reso- 
nant different-wavelengths propagation problem we are dealing with. 

III. TEMPORAL DYNAMICS OF LEVEL POPULATIONS AND RAMAN 
COHERENCE IN THE ADIABATIC APPROXIMATION (OPTICALLY THIN 

MEDIUM) 

In this section we study the temporal dynamics of populations and the atomic coherence 
in the given time- dependent field, assuming that the medium is optically thin. Based on 
that, Gi^2 will not depend on the coordinate z. One can show that condition (|l|) for Gaussian 
pulses reduces to G^^i ^ 1 when T2/T1 > With the first non-adiabatic correction, the 
solution of Eq. takes on the form: 

G2{r) Gi(r) 



G(r) ' ' " G(r) ' 
1 d{G2/G) 1 d{GJG) 



where G(r) = y^Gf (r) + Gi(r). 

The solutions for probability amphtudes are convenient to be represented as: 

ao = cos6'(r), 02 = — sin ^^(r), (7) 

where the mixing angle 0{t) is defined as tan^(r) = Gi{t) /G2{t) (we shall discuss its 
meaning later on). 

The expression for ai can be reduced to: 

ai = (G2G1 - GiG2)/G^ = e/G. (8) 



In adiabatic limit (0) |ai| = \0/G\ ^ 1 {0 = dO/dr), i.e. the population of the intermediate 
state |1) is close to zero all the time during the interaction with pulses. This also implies 
that the resonant absorption of the light pulses is weak (electromagnetically induced trans- 
parency) and the population is mainly distributed between the initial |0) and the final |2) 
states: 

|aoP + |a2p~l. (9) 

Equality (H) reflects the fact that atoms are trapped in the CPT state: Qcpt = {G2/G)ao — 
{Gi/G)a2 = ao cos 9 — a2 sin 6=1. This effect is responsible for the decrease in the resonant 
absorption of propagating pulses. Also the Raman coherence p2o = '^o'^2 occurs: 

P20 = -^sin(20) or pao = (10) 

Obviously, the maximum coherence (in absolute value) |p2o| = 1/2 is reached when 
9 = tt/4: {G^ = G^)- Figure 2 shows the temporal behavior of level populations |ao,2(T)P, 
the atomic Raman coherence |p2o(''")| and the mixing angle 6{t) for Gaussian pulses in the 
optically thin medium. 

The above results can be interpreted in terms of the three-dimensional vector model 
where vector variables ) and G = {G2, Gi, 0) (the torque vector) are introduced. 

Using these variables, we can rewrite (^) as 

a = G X a, (11) 

where the sign 'x' means the vector product. The solution of equation ([TT|) is vector a = 
{G2/G, —Gi/G,6/G). Components of the vector a coincide with adiabatic solution (^. 
Figure 3 illustrates the dynamics of vectors G and a in the three-dimensional vector model. 
The torque vector G moves in the 61-62 plane. And vector a having a small angle with 
respect to vector G {\9/G\ <^ 1) follows it. Thus one can see an absolute analogy with the 



adiabatic following in the case of a light pulse interacting with a two- level atom ||2^ . Such a 
simple picture can be observed only in optically thin media. In optically dense media, Gi^2 
and hence 6 become dependent on the z coordinate. 
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IV. SPATIAL EVOLUTION OF INTERACTING PULSES IN OPTICALLY 

DENSE MEDIA 



A. General case: unequal oscillator strengths {Ki / K2) 

The condition Ki 7^ K2 means that the probabihty of the |0) — |1) transition is not equal 
to that of the |2) — |1) transition. We note that in the ideal adiabatic limit ai = 0, and the 
pulses would not change their shape as they propagate in a medium that is optically thick 
for each of the pulses (see (^). However this is not the case. The non-adiabatic correction 
has to be introduced for real situations, which results in the induced dipole moments at the 
transitions |1) — |0) and |1) — |2) , and in the change of both pulses traveling in the medium. 
In order to attribute this effect to propagation of the interacting pulses in an optically dense 
medium, it is necessary to solve Eqs. (§) and @ in a self-consistent way. 

Use and (H) to rewrite field equations (H) in the following form 

From (p!2D, one can show that 

= K^Glir, z = 0) + K.Glir, z = 0), (13) 

i.e. K2G1{t, z) + KiG2{t, z) does not depend on the z coordinate. 

Equation (|T3|) describes the Manley-Raw relation, i.e. the law of conservation of the 
total energy density during propagation under CPT conditions. 

Using the definitions of 6 and G, we obtain the following expressions for 6*1,2 (t, -2) 

Gi{t,z) =G{t,z) sin 9{t,z), 

G2{t,z) =G{t,z) cos e{T,z). (14) 



Substitution (14) into the Manley-Raw relation gives the following expression for G{t, z): 
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r^2. . _ ^2,^^ sm^Ooir)) + K, cos^(go(r)) 

where Go('^) = + ^K'^'^O) ^'^'i ^o(t) = ^^o('r, 0) are the functions at the medium 

entrance 2 = 0. 

So the dynamics of level populations as well as the atomic Raman coherence and the 
evolution of the pulse shape are completely determined by the function 6'(r, z) that depends 
on both the time r and the coordinate z. Differentiating tan^ = G1/G2 on the z and using 
(|1^) we obtain the following equation for 6'(r, z): 



dl rn^dl^ (16) 
dr^ K{e) dz ' ^ ^ 

where K{e{T, z)) = Ki cos^ (^(r, z)) + K2 sin^ (^(r, z)). 

Eq. (|T6|) is similar to the equations describing nonlinear waves with the sharpening of 

the wave front during propagation [^]. The parameter u = G^/K can be treated as the 

"nonlinear" velocity. The nonlinear velocity can be described as u{t, z) = A{t)/ K'^{6{t, z)) 

where the first factor A{t) = GI{t)[K2 sin^ (^o(t)) + Ki cos^ (^o(t))] is independent of the z 

coordinate, and the second factor K{6{t, z)) is maintained along a characteristic of Eq.(|l6D 

6{t, z) = const. This allows to write down the characteristic curve equation in an evident 

form: 

1 



A{r')dr' (17) 



TO 



Here tq is the time at which the characteristic curve goes out of the medium boundary. 
The solution for 6{t, z) has the form: 

^(r,2) = 0o(ro,O) (18) 

Here tq is to be determined from Eq.(|l] 



B. The case of equal oscillator strengths {Ki = K2) 



In the case of equal oscillator strengths Ki = K2 = K the G(-2, r) function is not subject 
to changes during propagation ((^(r, 2;) = Gq{t) - the Manley-Raw relation for this case). 
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Therefore Eq. ([T6|) substantially simplifies: 

The solution of Eq. ( [T9|) can be written in the following form: 

e{T,z) = eo{z-\z{T)~z),o), (20) 

where Z{t) = JZoodr' G'^{0,t'), Z^^{t) is the inverse function of Z{t). 



It is not difficult to show that in this case G = \JG\{t, z) + G^t, z) coincides with the 
definition of dressed field pulses |]T2|,|T3[: G- = 00^2 — 02^1. Thus the pulses in our case 
can be identified as dressed field pulses (only at Ki = K2). It would be interesting to note 
that other combination G+ = a^Gi + 02^2 = (see also |12|). The concept of dressed field 
pulses cannot be applied to the case of Ki 7^ K2, but in both cases the Manley-Raw relation 
remains valid. 



C. The adiabaticity criterion 

The above results were obtained on the assumption that the adiabaticity criterion (|l]) 
(or \6/G\ ^ 1) remains valid during propagation of pulses. However it is not necessary the 
case. Therefore we investigate the adiabaticity condition with the account of propagation. 
Differentiating (0) on r, we can write the following expression for the adiabaticity criterion: 

e{T,z) _ deoG{T,z) K{e{T,z)) 

X [1 + sin (2^^o(ro))^^^2^^]-^ « 1 (21) 



As follows from (plD , the adiabaticity condition is destroyed {6/G ^ 00) when the factor in 
the square brackets tends to zero. Since sin (2^^o(to)) > in the entire range of change of 9 
(0 < < 7r/4), relation (^I]) is not fulfilled under the following conditions 

2{K^~K2)zd9o 



G'o(^o) 



sin(2^^o(ro)) 



{K,-K,)^>0, (22) 
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Evidently, condition (p2|) is not fulfilled at Ki = K2, and the adiabaticity criterion holds 
throughout the propagation process in this case. It should be mentioned that numerical 
analysis of Eq. (^) reveals that, generally speaking, there is a range of change of the 
q = K1/K2 parameter when the first condition in (|22| ) is not fulfilled either: g^m < q < qmax, 
qmin < 1, Qmax > 1- The valucs qmax and qmin depend on the ratio a = T2/T1 and the shape 
of pulses. For example, the q parameter for Gaussian pulses with T2/T1 = 3 may vary within 
the limit 0.40 < g < 1.35, as shown by the numerical simulation. 

In this case for the thick medium the adiabaticity criterion has the form: 

« 1 (23) 

It can be readily shown that condition ( p3D for Gaussian pulses is fulfilled when T2/T1 > -\/3 
and G2T1 ^ 1 (compare with the case of the thin medium). 

For Gaussian pulses (T2 > Ti) and the initial conditions ao(— 00) = 1, a2(— cxd) = 0, we 
have Oq > for —00 < tq < 0. With Ki > K2 {q > qmax) the adiabaticity criterion begins 
to break down at the leading edge: the mixing angle front becomes steeper. Condition ( pi]) 
gets destroyed at the trailing edge (where 6*0 < 0) when Ki < K2 {q < qmin)- Lets introduce 
a critical length Zc at which the adiabaticity condition (^) cannot be satisfied. Using (P^, 
one can obtain the following simple estimation for Zc for the case Ki 7^ K2: 

Let us now consider the experimental parameters: = lO^^cm"'^; 1/Ai^2 = 

10000, 20000(cm~^); 1/71,2 = 10,50ns are the relaxation times of levels |1) and |2), re- 
spectively; Ti = 0.1ns; G'?2^i = 20; a = T2/T1 = 3; fgio = 0.1 is the oscillator strength 
of the probe transition. From ( p^ ) we obtain Zc = ^O^Zq ^ 2cm, where zq is the linear 
absorption length. This estimation agrees with the results presented in Figs. 5-7. 

Using the characteristic equation (O) it is easy to find the area where the adiabaticity 



criterion breaks down. Characteristic curves for various values of q = K1/K2 are shown 
in Fig 4. The thickening of the characteristic curves means the sharpening of the mixing 
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angle front at a certain medium depth. At the point of intersection (see Fig. 4 b,c), 9 — > cxd, 
condition (|2T|) collapses. However there are no such points in the case when g^m < q < qmax 
(see Fig. 4a), i.e. here the adiabaticity criterion (pi]) is maintained during propagation of 
pulses. 

The above discussion is illustrated by Fig. 5 where the temporal behavior of the mixing 
angle 6 is presented for q = K1/K2 = 1,0.25,4 and for different normalized propagation 
lengths. On the Fig.5a-b the analytical results are presented (formulae (p!8| ) and (pO])). 
Figure 5c shows the numerical solution for 6 at Ki > K2. 

One can see that the evolution of the 6 parameter at Ki > K2 is different from that at 
Ki < K2- At Ki > K2 the adiabaticity condition fails for all values of z, beginning from 
the critical lengths Zc defined by Eq. (p^ . Here, the analytical theory does not apply at 
the very late stage of nonlinear wave propagation. In the case of Ki < K2 nonadiabaticity 
develops at the trailing edge (the front becomes steeper at a certain propagation length), 
but it does not go deeper into the medium. A good agreement between the analytical and 
numerical solutions for ^ at i^i < K2 is observed over the entire propagation length. This 
leads us to conclude that the interaction adiabaticity is fairly sensitive to the ratio of the 
oscillator strengths of the transitions interacting with the pulses. 

D. Discussion of results 

The solutions obtained have been used to analyze the temporal dynamics and spatial 
behavior of the propagating EIT pulses and the atomic coherence for various oscillator 
strength ratios. Figure 6 illustrates the temporal and spatial evolution of normalized Rabi 
frequencies of the both pulses 5'i,2(t) = ^1^2/ \J + as they propagate inside the 
optically thick medium. The temporal evolution of pulses can be seen to depend on the ratio 
between the transition oscillator strengths. In the case of Ki < K2 both pulses undergo 
reshaping as they propagate in the medium (Fig. 6a and b). The probe pulse is gradually 
depleted and the coupling gets stronger. Note that the pulse shape at the initial stage of 
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propagation shows very little change along the length of the medium, which may exceed 
the linear absorption length. Complete reemitting of the probe pulse into the coupling one 
during propagation is possible. Using the Eq. ( p]7D one can obtain the following expression 
for the maximal distance Zm on which the probe pulse propagates into the medium: 

1 



K\e,{-oc)) J-oo 

For the Gaussian pulses we have 



A{T)dT. (25) 



Here a = T2/T1. 

For the same parameters as in the previous sub-section, we have Zm = 4.8 ■ lO^^o ~ 9cm 
for K1/K2 = 1 and Zm = 8.4 ■ lO^^o ~ 16cm for K1/K2 = 1/4. These values agree with the 
results shown in Figs. 5-7. 

An interesting feature of the spatial distribution of the probe pulse is illustrated in Fig. 
6a,b: in some areas of the medium, the field in the tail of the pulse is different from zero. 
This is believed to be the result of the spatial compression of the probe caused by the slowing 
down of the group velocity of the probe pulse EH . 



In the case of Ki > K2, the adiabaticity condition is maintained over the Zc range, which 
can also be much longer than the length of the linear absorption. The leading edge of the 
probe pulse undergoes gradual depletion and the pulse amplitudes display only small changes 
in that range. Outside that range, the pulse splits into several peaks. This occurs due to 
the nonadiabaticity of interaction (see Fig. 6c where numerical simulations are presented 
for (71,2 (t)- Here the first three curves correspond to the adiabatic interaction). 

The pulse evolution, described above, is determined by spatial and temporal behaviour 
of the atomic Raman coherence P20 shown in Fig. 7 (see also |0). Figures 7a,b reveal 
an unusual spatial and temporal behavior of the atomic coherence, which we interpret as 
the slowing down, stopping and localization of the atomic coherence in the medium. The 
probe pulse is transferred into and stored in the collective atomic excitation under the 
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control of the coupling. One can say that a phase grating is created in the atomic medium. 
The phase grating is preserved throughout the entire period of relaxation of the atomic 
coherence. During that period of time, the information stored in the atomic excitations 
can be transferred back to the radiation using another coupling pulse of the same or of a 
different frequency. 

We find these effects to be similar to the ones predicted and demonstrated in P, p!o| , ^ 



but for some differences as indicated below. In our case, both the coupling pulse and the 
probe pulse are strong. Unlike |^ , we used the boundary conditions for the pulse envelopes 
yielding the time distribution of pulses at the medium boundary, z = 0. We believe, this 
condition is more natural than the one used in |^ where the authors use the probe pulse 
distribution in the medium at fixed time as the initial condition. Also they do not take into 
consideration evolution of the coupling pulse. The effect takes place in the case Ki < K2 
and is not observed in the case of Ki> K2. 



V. CONCLUSION 

The propagation of two short overlapping pulses with durations T2 > Ti in optically thick 
three-level media under CPT conditions has been studied for the general case of unequal 
transition oscillator strengths. An analytical solution has been obtained for the set of reduced 
wave equations under the adiabatic following condition. Also it has been shown how spatial 
evolution of pulses depends on the oscillator strength ratio. 

The condition of adiabaticity provided at the medium entrance preserves for any value of 
propagation lengths if Ki < K2 {q < Qmax) and breaks down at Ki > K2 {q > Qmax)- In the 
range q < qmax, the probe pulse is completely depleted and reemitted into the coupling pulse 
during propagation. This is not possible in the case of g > qmax- It has been established 
that to provide for the adiabaticity condition in an optically thin medium, the restriction 
T2/T1 > \f2 has to be ensured, whereas in a thick medium T2/T1 > ^/S. 

We also have studied the spatial behavior of the atomic coherence P20, which plays a 
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significant role, for example, in nonlinear mixing processes. It has been found that a strong 
coherence can be maintained over a length equal to several hundreds of thousand of one- 
photon absorption lengths during propagation. The effect of localization of the atomic 
coherence is demonstrated. 
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FIGURES 

FIG. 1. (a) The three-level Lambda-type system coupled by two resonant pulses with 
Rabi frequencies Gi and G2. (b) The shapes of probe Gi and coupling G2 pulses at the 
medium entrance. 

FIG. 2. The time evolution of (a) level populations |ao,2(''")P and the atomic Raman 
coherence |p2o(t)|, (b) the mixing angle 6{t) in an optically thin medium for the Gaussian 
pulses. G?Ti = G^Ti = 20, (a) T2/T1 = 3; (b) T2/T1 = ^2 (A), ^3 (B), 3 (C), 10 (D). 

FIG. 3. The vector model of adiabatic interaction of two short pulses with a three-level 
Lambda-type system. 

FIG. 4. The characteristic curves for the Eq.(|16D: (a) K1/K2 = 1, (b) K1/K2 = 0.25, 
(c) Ki/K2 = 4. 

FIG. 5. The time evolution of the parameter 6 for different relationships between Ki 
and K2 at different propagation lengths, (a) Ki = K2, (b) K1/K2 = 0.25, (c) K1/K2 = 4. 
G°2^i = 20, T2/T1 = 3. Here and in all the other figures the time r is measured in the 
units of the pulse duration Ti, and the propagation length z of pulses in the medium is 
measured in the units of the length of linear absorption of the probe radiation determined 
in accordance with the Beer's low. In (c) the numerical solution for the case Ki/ K2 = 4: is 
presented. 

FIG. 6. The time evolution of the normalized Rabi frequencies (71 2 = ^1^2/ \/g^~+~G^ 
of the probe and coupling pulses for different relationships between Ki and K2 at different 
propagation lengths, (a) Ki = K2, (b) K1/K2 = 0.25, (c) K1/K2 = 4. = 20, T2/T1 = 3. 

FIG. 7. The time evolution of the atomic Raman coherence for different relationships 
between Ki and K2 at different propagation lengths, (a) Ki = K2, (b) K1/K2 = 0.25, (c) 
K^/K2 = 4. G? 2T1 = 20, T2/T1 = 3. 
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